Precipitation and river discharge are critical components that define hydrological regimes of watersheds. Recently, this region has experienced drought conditions and water resource managers need detailed information about changing flow conditions to better plan for decreased streamflow. The purpose of this report is to provide an analysis of the long term trends in precipitation and discharge in the Breitenbush River watershed in the Oregon Cascades, USA. In this analysis, daily streamflow data will be used to detect changes in the streamflow regime by calculating the annual daily mean for multiple periods from 1932 to 20202. The change in mean will be expressed as percent, slope (%/year), and standard deviation of discharge. R tools from the US Geological Survey will be used to import data and perform the long-term analysis. These tools include dataRetrivial and EGRET.
Declining precipitation as snowfall, is causing runoff to occur earlier and more rapidly (Feng and Hu, 2007). This can lead to larger peakflow due to flashy rainfall events. In addition, low flow conditions can be effected by changes in precipitation phase change because snowpack releases water gradually into summer. This can can impact the water resiliency of communities and ecosystems that rely on late season base flow (Rumsey et al., 2015). The Oregon Cascades are a relatively low elevation mountain range and is vulnerable to changes in surface temperature and snow pack reduction. The Breitenbush River subwatershed (HUC 17090005) is located in the North Santiam Watershed in Marion County, Oregon. The climate is warm-summer mediterranean and receives about 59 inches of rain annually. The subwatershed also receives about 146 inches of snowfall annually. Most of the precipitation falls during the winter months. Streamflow monitoring has been ongoing since 1932 and a USGS stream gaging station (44.75262), -122.129) is the primary data source for this study. The upstream drainage area is 279.7 square miles.
Figure 1. Location of Breitenbush subwatershed with stream gage location and upstream drainage area.
Packages developed by the U.S. Geological Survey were used to determine long-term trends in streamflow. The first is dataRetrieval package (De Cicco et al., 2018) which aids in the reproducability of the results. This package enables rapid acquisition of hydrological data from the National Water Information System which is database containing historical and instantaneous time-series. The entire daily discharge record was retrieved; the use of daily discharge values for computation of historical streamflow statistics is common practice in hydrology. Due to a large data gap from 1987 to 1998 when the gaging station was not active, as a result two periods were analyzed using the EGRET package (Hirsch & Cicco, 2015). This package allows for exploratory data analysis of trends in river data, including water quality trends.
Two aspects of stream flow trends were investigated. The first is the occurrence and variability of 7-day low flow which is the average discharge during the 7 days with the lowest flow for the period specified. Low flow trends are correlated with anthroprogenic water use and precipitation variability (Hammon and Flemming, 2021). Second, maximum daily trends were analyzed during the winter months (November through February).These streamflow metrics were used to assess streamflow response to changes in the precipitation regime.
Maximum daily flows increased substantially through from the beginning of the streamflow record until the late 1960’s when this trend flat lined and eventually began a gradual negative trend that continues to the present day. A similar, but more defined trend is present when 7-day minimum flow is investigated. In this case, an upward trend of 24 percent is observed until the late 1960’s when this trend reverses. A negative 5.5 percent decline in 7-day minimum flow is observed from about 1968 to 1987 and a negative 15 percent decline is record between 2000 and 2019.
Figure 2. Plot of maximum daily streamflow metrics for 1936 to 1987.
Figure 3. Plot of 7-day minimum flow from 1934 to 1987.
BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
Season Consisting of Jul Aug Sep
maximum day
Streamflow Trends
time span change slope change slope
cms cms/yr % %/yr
1932 to 1950 1.5 0.085 14 0.8
1932 to 1968 2.6 0.072 24 0.68
1932 to 1987 1.9 0.034 17 0.32
1950 to 1968 1.1 0.059 8.8 0.49
1950 to 1987 0.33 0.009 2.7 0.074
1968 to 1987 -0.74 -0.039 -5.5 -0.29
Table 1. Trend analysis for maximum daily streamflow, 1932 to 1987.
BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
Season Consisting of Jul Aug Sep
7-day minimum
Streamflow Trends
time span change slope change slope
cms cms/yr % %/yr
1932 to 1950 0.43 0.024 14 0.76
1932 to 1968 0.59 0.016 19 0.52
1932 to 1987 0.37 0.0068 12 0.21
1950 to 1968 0.16 0.0088 4.4 0.24
1950 to 1987 -0.06 -0.0016 -1.7 -0.045
1968 to 1987 -0.22 -0.012 -5.8 -0.31
Table 2. Trend analysis for 7-day minimum streamflow, 1932 to 1987.
Figure 4. Plot of maximum daily streamflow, 2000 to 2019.
Figure 5. Plot of 7-day minimum streamflow, 2000 to 2020.
BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
Season Consisting of Nov Dec Jan Feb
maximum day
Streamflow Trends
time span change slope change slope
cms cms/yr % %/yr
2000 to 2010 1.1 0.11 1 0.1
2000 to 2019 -1.8 -0.093 -1.6 -0.082
2010 to 2019 -2.9 -0.32 -2.6 -0.28
Table 3. Maximum daily streamflow, 2000 to 2019.
BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
Season Consisting of Jul Aug Sep
7-day minimum
Streamflow Trends
time span change slope change slope
cms cms/yr % %/yr
2000 to 2010 -0.28 -0.028 -7.9 -0.79
2000 to 2019 -0.54 -0.028 -15 -0.8
2010 to 2019 -0.26 -0.029 -7.9 -0.87
Table 4. 7-day minimum streamflow, 2000 to 2019.
The Breitenbush River has seen periods of increasing streamflow metrics during the first half of the 20th century, but recent trends indicate a reversal of these trends. Water resource managers downstream of this basin should be aware of the severity of this negative trend, because it may have implications for water security. Communities that rely on summer time low flows for irrigation could faces reduced water resource resiliency. As Oregon enters another year of record drought, it is important to note the historical trends do not indicate that we will return to normal conditions in the near future.
De Cicco, L. A., Lorenz, D., Hirsch, R. M., & Watkins, W. (2018). dataRetrieval: R packages for discovering and retrieving water data available from US federal hydrologic web services. Reston, VA, DOI, 10, P9X4L3GE.
Feng, S., & Hu, Q. (2007). Changes in winter snowfall/precipitation ratio in the contiguous United States. Journal of Geophysical Research: Atmospheres, 112(D15).
Hirsch, R. M., & De Cicco, L. A. (2015). User guide to Exploration and Graphics for RivEr Trends (EGRET) and dataRetrieval: R packages for hydrologic data (No. 4-A10). US Geological Survey.
Rumsey, C. A., Miller, M. P., Susong, D. D., Tillman, F. D., & Anning, D. W. (2015). Regional scale estimates of baseflow and factors influencing baseflow in the Upper Colorado River Basin. Journal of Hydrology: Regional Studies, 4, 91-107.